The 130 GeV gamma-ray line and generic dark matter model building constraints 
from continuum gamma rays, radio and antiproton data 
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An analysis of the Fermi gamma ray space telescope data has recently revealed a resolved gamma- 
ray feature close to the galactic center which is consistent with monochromatic photons at an energy 
of about 130 GeV. If interpreted in terms of dark matter (DM) annihilating into 77 (jZ, 7/1), this 
would correspond to a DM particle mass of roughly 130 GeV (145 GeV, 155 GeV). The rate for these 
loop-suppressed processes, however, is larger than typically expected for thermally produced DM. 
Correspondingly, one would generically expect even larger tree level production rates of standard 
model fermions or gauge bosons. Here, we quantify this expectation in a rather model-independent 
way by relating the tree level and loop amplitudes with the help of the optical theorem. As an 
application, we consider bounds from continuum gamma rays, radio and antiproton data on the 
tree level amplitudes and translate them into constraints on the loop amplitudes. We find that, 
independently of the DM production mechanism, any DM model aiming at explaining the line 
signal in terms of charged standard model particles running in the loop is in rather strong tension 
with at least one of these constraints, with the exception of loops dominated by top quarks. We 
stress that attempts to explain the 130 GeV feature with internal bremsstrahlung do not suffer from 
such difficulties. 

PACS numbers: 



I. INTRODUCTION 

A very promising way to search for dark matter (DM) 
particles is to look for the pronounced spectral signatures 
in cosmic gamma rays which are expected from DM anni- 
hilation or decay [l|. A (quasi-) monochromatic gamma- 
ray line, in particular, results if DM annihilates or decays 
into final states including a photon (such as 77, 7Z or 7/1) 
Q or, in particular for Majorana DM particles, due to 
internal bremsstrahlung (IB) photons for ffj final states 
Such a feature is usually considered a smoking gun 
signature for particle DM as there are no known astro- 
physical processes that produce mono-energetic gamma 
rays at high energies. 

The recent discovery 0, Q of a tentative line-like spec- 
tral feature in the data of the Fermi Gamma-Ray Space 
Telescope @ (independently confirmed in Ref. Q) has 
therefore prompted considerable interest in the field, in 
particular after a sophisticated spatial template analy- 
sis by Su & Finkbeiner has confirmed its existence with 
an impressive global significance of slightly more than 5 a 
Q . Instrumental systematics seem unlikely as a cause for 
this gamma-ray line at an energy of around 130 GeV @ 
and no compelling alternative explanation to DM annihi- 
lation has been brought forward so far (but see Ref. [l(| ) ■ 
However, there also exist a few indications which even- 
tually may, if confirmed, disfavor the DM interpretation, 
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so more data is certainly needed to settle this issue. For 
more details about the current status of the signal and 
possible implications for DM models, including an exten- 
sive list of references, we refer the reader to a dedicated 
review [Til ]. 

One possible worry is that the required annihilation 
rate is considerably larger than typically expected for 
thermally produced DM [H, [l]| , which means that one 
may have to abandon this theoretically very appealing 
mechanism for DM production in the early universe. 
Even in that case, however, a large annihilation rate into 
77, 7Z or 7/1 final states usually implies a much larger 
annihilation rate of DM particles into pairs of standard 
model (SM) fermions or weak gauge bosons - simply be- 
cause the former processes can only happen at the 1- 
loop level, while the latter happen at tree level. It has 
been noted that these annihilation products could there- 
fore rather easily violate existing bounds on continuum 
gamma rays fi#fl7t or cosmic-ray antiprotons [ll| in typ- 
ical scenarios where DM consists of weakly interacting 
massive particles (WIMPs). In fact, it has already been 
claimed [l6| that the leading DM candidate, the lightest 
neutralino in supersymmetric theories [HI , is ruled out as 
an ex plan ation of the observed line signal (see, however, 

Refs. inm-si)- 

Here, we revisit these arguments in an as systematic 
and model-independent way as possible. As sketched in 
Fig. [TJ we use the optical theorem to relate the expected 
annihilation rates at tree-level to the imaginary part of 
the loop-level amplitude that is supposed to account for 
the line signal (extending, and in a certain sense turn- 
ing around, the recent analysis by Abazajian et al. |28j|). 
We then confront these trec-lcvel annihilation rates, on a 
channcl-by-channel basis, to constraints arising from the 
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Figure 1: Graphical illustration of Eq. i.e. the optical theorem applied to the case of DM annihilation into yX, where X = 
7, Z, h: all intermediate-state particles I that can be put on-shell contribute to the imaginary part of the loop-amplitude. This 
allows to relate the latter to the tree- level annihilation rate of DM into SM particles with msM <i r n x , where the corresponding 
ii-channel diagram (as well as the 4-point interaction for W + W~ — > ~yX) is not shown explicitly in the figure. For yZ and 
7/1 final states dominated by fermion loops, another possible way of putting the intermediate particles on-shell - compared to 
what is shown on the right hand side - results from 'cutting' the loop diagram along the dotted rather than the dashed lines. 



non-observation of excesses in continuum gamma rays 
from the galactic center and dwarf galaxies, as well as 
antiproton and radio data. 

Our main result is that the expected large tree-level 
annihilation rates are indeed both a generic and serious 
challenge to any DM model trying to explain the observed 
130 GeV feature by a loop that involves SM particles. 
Given that DM usually couples directly to SM particles 
in commonly adopted WIMP frameworks like supersym- 
metry, this has severe implications on concrete model 
building attempts. We classify the (few) possible loop- 
holes to the above statement and stress that IB (which 
upcoming experiments could distinguish sufficiently well 
from monochromatic photons [29| ) may well be a more 
likely explanation of the signal in terms of DM annihi- 
lation. Alternatively, our results seem to imply that the 
dominant coupling of DM to charged states must be ex- 
clusively to new particles which are heavier than the DM 
particles themselves - a possibility which may not be very 
common in typical WIMP scenarios (aiming simultane- 
ously at addressing the fine-tuning issues of the SM) but 
which has seen a lot of interest in response to the observed 
photon excess at 130 GeV (see, e.g., Refs. 0^27)]). 

This article is organized as follows. We start in Section 
|TT] by a short review about the optical theorem and how 
one can use the strength of the potentially observed line 
signal at 130 GeV to derive a lower bound on the expected 
DM annihilation rate at tree level. Constraints on tree- 
level annihilation rates from both continuum gamma-ray, 
antiproton and radio observations are then discussed in 
Section IIIIl We combine these results in Section IIVI to 
constrain the maximally allowed fractional contribution 
of the imaginary part of the loop-amplitude to the full 
annihilation cross section into 7 V (with X = 7, Z, h) 
and discuss implications for DM model building. In Sec- 
tion |Vl we present our conclusions and an outlook. In 
the Appendix, we collect some technical details about 
our calculations and the application of the optical theo- 
rem to the particular processes we are interested in here. 
Throughout the paper we use natural units, Co = K = 1, 
unless indicated explicitly. 



II. FROM LOOP TO TREE-LEVEL PROCESSES 

Our main goal is to establish a relation between the 1- 
loop level cross sections a 7 x for dark matter annihilation 
into 7V final states and the tree level cross sections aj 
to pairs of Standard Model or other particles which form 
intermediate states denoted by I for the production of 
this final state. To this end, we express the unitarity, 
S^S = 1, of the 5-matrix 5" = 1 + iT in terms of the 
transition matrix T, 

-i{T-T r ) = T^T. (1) 

Sandwiching this equation between an initial and a final 
state \i) and |/), respectively, yields 

S[(/|TH)] = i^</|Tt|X>(X|TH) , (2) 
x 

where the sum runs over all on-shell intermediate states 
I X) whose total four momenta equal the initial state four- 
momentum p, in our case consisting of two dark matter 
particles Xi an d whose total angular momentum equals 
the one of the initial and final state (J, M). Note that 
this assumes (/|T T |i) = (f\T\i)* - as is, e.g., satisfied in 
the angular momentum basis if the theory is invariant 
under time- reversal [28j | , which is the situation we are in- 
terested in here (the same would be true when using a 
helicity basis). Obviously, \X) must also share all other 
conserved quantum numbers with \i) and |/), such as in 
particular parity. We now assume that this sum only 
contains states consisting of two free particles of type /. 
Exploiting energy momentum conservation in the cen- 
ter of mass frame all relevant two particle states can be 
uniquely characterized by the total spin S of the two par- 
ticles, its projection ms on some given coordinate axis, 
and the unit vector p x of one of the particles. Thus for 
the angular part of the wave functions for such states we 
can write 
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\X) 



|pi, S,n%s}i = ^ i( L , m L,S,ms\pi,S,ms)i\L,m Ll S,m s )i = ^ Y L , mL (pi)\L, m L , S,ms)i, 



r 



(3) 



where we have expanded this state in terms of the 
cigenstatc \L,mi,, S,ms)i of orbital angular momen- 
tum and spin and we have used that the expan- 
sion coefficient is simply the spherical harmonic func- 



J 



tion, i(L,m Ll S,m,s\pi,S,m s )i = Y Ll m L {pi)- Ap- 
plying (/IT^I from the left to Eq. §3§ and inverting 
to (f\T'\L, rtiL, S, ms)i by using the orthonormality 
J d 2 piYZ imi (pi)Y L2 , m2 (pi) = S LuL2 S mum2 gives 



(f\T^\L,m L ,S, m s )i 
(J\T*\vi,S,ms)i 



J d 2 P iy iimi (pi)(/|Tt| Pl ,5, m s )j = (J, M\L, m L , S, m s )iM}^ LiS) (a) 

Y Lr, lL (Pl)( J i M \ L > m L,S,™s)lM}^ {I ^ s) (s). 



(4) 



L,m L 



r 



Here we have used the fact that if the transition ma- 
trix is a scalar the amplitude (f\T'\L, raj,, S, ms)i has 
to be proportional to the Clebsch-Gordan coefficient 
(J, M\L, m,L, S, mg) and we denoted by M.%-^ti t L,S){ s ) 
the remaining factor which is independent of the pro- 
jections mg and tul and can only depend on the total 
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angular momenta L and S and on the Lorcntz invari- 
ant s = p 2 = (pi+p-2) 2 - Eq. (0| provides the relation 
between the amplitude (f\T^\pi, S,ms)i that is calcu- 
lated directly from the Feynman rules and the amplitudes 
■M-i^f(i,L,s) ( s ) that are more convenient for the following 
analysis. For Eq. @ we can now write 



2-Ei (2tt) 3 J 2E 2 (2?r) 



I 



(/|Tt| Pl , ft m s ) II {p 1 ,S, ms\T\i)(2n)H\ Pl +pa - p) . (5) 



r 



Inserting the second equation in Eq. Q and its ana- 
logue for /(pi, ft ms\T\i) and performing the momen- 
tum integration in Eq. ((5]) the orbital part drops out 
due to the orthogonality of the spherical harmonic 
functions and the spin dependent part simply gives 
T, m \{J,M\L,m,S,M - m) | = 1 due to the unitarity 
of the Clebsch-Gordan coefficients. One ends up with 
the reduced velocity 



ft 



ft +ft 



(6) 



where ft = (p^/I?, = Pcm/^x are ^ ne velocities of the 
two intermediate state particles with pcm the momentum 
in the center of mass frame. This reduces to ft = ft = ft 
for equal mass particles. Overall, one gets 

9 [M^f(s)} -E#2E M^ {I , L , s) (s)M% {I . L . s] ( S ) , 

/ Z n L,S 

(7) 

where the sum over L and S runs over all combinations 
which can combine to J and are allowed by possible other 



conserved quantum numbers carried by the intermediate 
state, such as parity. 

As schematically shown in Fig. [TJ Eq. ([5]) relates two 
tree level amplitudes to a 1-loop process by cutting the 
loop along the dashed lines and putting the intermediate 
states on-shell. Note that for the special case \ i) = \f) the 
right hand side of Eq. is essentially the total cross sec- 
tion into all possible intermediate states which 
is thus related to the imaginary part of the forward scat- 
tering amplitude (i\T\i) = .Mi-n(p). This is known as 
the optical theorem. 

The transition rates from a given initial state 

\i) to any state / consisting of two particles with not 
necessarily equal masses are proportional to the phase 
space integral on the right hand side of Eq. ([5]) taken for 
|/) = With Eq. (J7J, and including the appropriate 
normalization factors, this reads 



i->7 



(s) (X va i ^ I (s) = 



2 5 7T 2 



L,S 



\M lMI . L . s) {s)\ 2 . (8) 
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While wc keep here all expressions fully general, we will itics and therefore could in any case not account for the 



later mostly consider the v = 2/3j — > limit and thus 
restrict ourself to s-wave annihilation, motivated by the 
fact that p-wave contributions are strongly suppressed 
with a factor of v 2 ~ 10~ 6 for typical galactic DM veloc- 



observed line signal. Squaring Eq. ([7]), using the Cauchy- 

Schwarz inequality |E; a A| 2 < (J2i \ a i\ 2 ) (J2i M 2 ) and 
combining with Eq. ([8]) yields 



J 



i-s-7 



< 



N f (3f 



V 



2 12 7T 4 ^ 
I.L.S 



Pi\M { 



I,L,S)- 



2 6 7T 2 



(9) 



r 



where 



\Mi- 



< 1 



(10) 



The symmetry factor Nf introduced above is 1/2 for two 
photon final states and 1 otherwise; note that such a 
factor does not appear in Eq. © because we are only 
interested in charged intermediate particles. 

If there is only one intermediate state (/, L, S) or if 
one of them dominates, Eq. @ turns into an equality. 
We note that if the absolute value of the amplitudes 
■^i~>(i,L,s) an d M. j_>(7 l s) d° n °t correlate, the limit 
Eq. @ is very conservative and one should resort to the 
original sum in Eq. (JT)) to derive limits. In order to pro- 
ceed in that case, one has to specify an (effective) La- 
grangian which describes the interaction between \ and 
J, see Appendix IA 21 

Eq. ([9]) is our master equation. We will apply it to the 
case where the final state |/) is 77, 7Z or 7/1, summed 
over all polarization states with the same J p value as 
the initial state i of the two annihilating DM particles, 
and the intermediate states \I) consist of a particle-anti- 
particle pair. Since the angular momentum orientation 
independent part M.(i^L,s)^f of the tree level amplitude 
defined in Eq. ((4]) can straightforwardly be calculated 
from the Standard Model or its extensions and the tree 
level cross section is constrained from observations 
to be discussed in Sect. IIII1 we will get constraints on 
r"i_>/ times the 1-loop cross section cr^f giving rise to 
gamma-ray lines; for this quantity, we will sometimes use 
the short-hand notation. 



(11) 



where the right-hand side is understood to be summed 
over all polarization states as explained above. 

For a scalar initial state, J = 0, the sum reduces to a 
sum over L = S only; this is for example the case for Ma- 
jorana DM in the v — > limit we are interested in here. 
In this case the initial state has also negative parity, im- 
plying that fermionic (bosonic) states / must combine 
to L = S = 2n (L = S = 2n + 1), with n = 0,1,2.... 
Table |T] shows the corresponding ratios of Eq. © for all 
relevant channels. The results for scalar DM, where the 



initial state has positive parity, are given in Tab.llll here, 
L = S = 2n (2n + 1) for bosonic (fermionic) states / - 
which implies that there are two independent degrees of 
freedom for W boson loops (see also the discussion in 
Appendix IA In the following, we will limit our dis- 
cussion to the Majorana case but note that an extension 
of our analysis to DM particles with any spin properties 
in the initial state is straight-forward with the formalism 
presented here [30}. For more details about the calcula- 
tion leading to the results presented in Table HI we refer 
the reader to Appendix [5] We checked explicitly that 
the two methods described there give consistent results, 
and that the results for 77 final states agree with those 
reported in Ref. [H]. 

Let us finally point out that in the case of "/Z and 
7/1 final states dominated by fermion loops, there is a 
contribution to the imaginary part of the cross section, 
at the same order in the coupling constants, that is not 
captured in Eq. §§§ above. This is related to the fact 
that both Z and h are unstable, such that there is an- 
other way of cutting through the loop diagram shown 
in Fig. [T] (as indicated by the dotted line). While the 
possibility of DM annihilation to such three-body (ff"f) 
final states is obviously not taken into account in our for- 
malism leading to Eq. ([9]), we note that those processes 
are 0(a em )-suppressed with respect to the tree-level an- 
nihilation (to //) and therefore only marginally affect 
the general limits on the annihilation rate that we will 
discuss in the next Section. While the contribution to 
the imaginary part of the loop amplitude still may be 
sizable, explicit calculation for Majorana DM along the 
lines explained in Appendix IA 21 shows that those contri- 
butions always come with opposite signs as compared to 
the contributions from Eq. ©. Neglecting them will thus 
lead to an over- estimation of the imaginary part and thus 
conservative bounds on the quantity defined in Eq. pip . 



III. CONSTRAINTS ON THE ANNIHILATION 
RATE 

At tree-level, DM annihilates or decays into SM par- 
ticles which then (except for electrons and neutrinos) 
fragment and/or decay into stable particles like pho- 
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IVlu 1 ' ' I CL1.1.CL 

DM 


'final' state T 


(<™)tree 


(cru)tree 


(CTu)trcc 


XX -> // 


(0,0) 


2^-^-^ [tanh /3/J 


c OS ^ w ^ [ 2 Qsm 6 W \ ^[tanh /3jJ 





xx -> w+w~ 


(1,1) 


2aL/3/ [tanh- 1 /?/] 2 


4a cm acos 2 Owfiffii [tanh^ 1 /?/] 2 






Table I: For Majorana DM, the initial state has J = 0~ in the zero- velocity limit, leaving only the above stated spin and orbital 
angular momentum combinations for charged 'final' state particles /. The quoted relations between tree-level and loop-level 
annihilation rates, with /3f = 1 — m%/4m x and 0j — 1 — m 2 ^ w ^/m x , then follow as described in the text, c.f. Eqs. (JUHTTJ), and 
Appendix 1X1 



Scalar 
DM 


'final' state / 
(L,S) 


3(o"«)t-, 




3(<™) T // 
(CT-u)tree 


xx // 


(1,1) 


ZSr-t-Pi [tanh 1 /9i] a 


5-2 a— a T o — Q sm "W 

cos^ 6\yPfPjm^ 12 ^ J 

x^tanh-^-^lfi-L^^!)] 2 





xx -> w+w~ 


|t)=^/f|0,0) 
+ 7fl 2 > 2 > 

io=\/ii 2 ' 2 ) 
-^i°.°> 


2«L(l-^|)^[tanh-V/] 2 
aL^^Itanh- 1 /?/] 2 


cos 2 0W [(^2 + ^ tanh -l ^ 






2/3/0/ m 


« L V" 4m^, 8mi, 7 ^ J 
'n «2 m| W tanh-^j , \ I 2 
2 ft 2^,J(, ^ !JJ 



Table II: Same as Tab. U but for Scalar DM (where the initial state has J = + in the zero- velocity limit). Note that in this 
case there are two independent degrees of freedom for the W-boson 'final' states /; here, we state the results for longitudinal 
and transverse states (see also Appendix I A II for more details). 



tons (mostly via 7r° — > 77), antiprotons and electrons 
or positrons. In this Section, we will discuss the var- 
ious constraints on the decay or annihilation rate that 
derive from the non-observation of excesses in gamma- 
rays, cosmic ray antiprotons and radio data. We will 
focus on DM masses that could explain the observed line 
feature at E^ 30 = 129.8^ 3 4 2 GeV Q, which in the case 
of annihilating DM implies 



£130 



(12) 



where X = y,Z,h denote the most interesting possibil- 
ities (thou gh in principle X may also be a new neu- 
tral state [31|). Decaying DM would require m x = 



£130 

7 



2m 2 x , but we will not consider this 



option in the following because it is not supported by the 
observed angular distribution of the signal [IJ (3 ■ 

For all of this section, we will for consistency use an 
Einasto profile for the DM distribution in the Milky Way, 



2 r a 

Px( r ) = Poexp ( -- — 



(13) 



with parameters chosen such as in Ref. 0, i.e. a = 0.17, 
r s = 20kpc and po = 1.05 • 10 4 GeVcm 3 (which results 
in p x (Rq) = 0.4 GeV cm 3 at the position of the sun, 
Rq = 8.5 kpc). Such a profile is not only favored by 
A^-body simulations of gravitational clustering [32[ but 
also provides a very nice fit to the line observation (ll| . 



A. Continuum gamma rays 

Dwarf spheroidal galaxies exhibit mass-to-light ratios 
of up to 0(1000), which are the largest values observed 
in astronomical objects. If located at not too large dis- 
tances from the sun, this makes them ideal targets for 
DM detection because the gamma-ray emission expected 
from astrophysical processes is negligible [H, Hj]. In 
fact, the non-observation of any gamma-ray signal from 
dwarf galaxy satellites of the Milky Way by Fermi places 
the currently strongest constraints on secondary photons 
from DM annihilation [35[. We collect these constraints 
for all relevant channels in Tab. lIIIi conservatively choos- 
ing the weakest bound in the ranges for the DM masses of 
interest, as defined by Eq. (|T2|) . We note that constraints 
on the channel are not explicitly provided in Ref. [35| . 
However, the photon spectrum from is very similar to 
the one from /j, , which means that a good estimate for 
bounds on e final state is obtained by rescaling the /i 
bounds by the number of photons produced per annhila- 
tion, N 7 = J^ydEjdN^/dE^, which in both cases is 
given by final state radiation [36[. 

The Fermi dwarf limits assume a DM density profile 
that scales as r _1 in the innermost part. Such an NFW 
profile is also consistent with numerical A^-body simu- 
lations [37} and only slightly steeper in the central part 
than our reference profile of Eq. (p~3|) . While kinemat- 
ical data presently cannot distinguish between shallow 
or (possibly very) cuspy central profiles in dwarf galax- 
ies, the resulting limits on the DM annihilation rate can 
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change sizably. For that reason, a more conservative ap- 
proach relies on considering instead limits on the annihi- 
lation rate that derive from continuum gamma rays in the 
galactic center region itself; since we will eventually be 
interested in the ratio of photon fluxes for the continuum 
and line component, this also minimizes the astrophysical 
uncertainties related to the DM density profile. 

In Tab. IIIII we therefore also present a compilation of 
the corresponding limits from a recent analysis of a re- 
gion with \b\ < 5° and \l\ < 5° around the galactic center 
[17| . Those limits were derived from a slightly different 
DM profile than in Eq. (fT3"|); we therefore divided them 
by a factor of 1.4, taking into account that they should 
roughly scale with the quantity J = J An dfl J t ds p 2 (r) 
(where Af2 corresponds to the angular region considered 
and l.o.s refers to a line-of-sight integration). Note that 
the limits on light lepton final states are actually stronger 
than for the dwarf analysis because of inverse Comp- 
ton scattering of high-energy e ± , which upscatters low- 
energy photons of the rather dense interstellar radiation 
field to gamma-ray energies. 



B. Antiprotons 

The propagation of antiprotons can very well be de- 
scribed in simple phcnomcnological diffusion models, the 
parameters of which are strongly constrained by other 
cosmic ray data like in particular the boron over car- 
bon ratio B/C [38| . In this way one can predict the 
astrophysical background (dominated by secondary an- 
tiprotons from cosmic ray proton collisions with interstel- 
lar hydrogen) with remarkably small uncertainties. The 
flux of primary antiprotons from DM annihilation, on 
the other hand, is subject to greater uncertainties; these 
are mostly related to the vertical size L of the magnetic 
halo which is directly proportional to the annihilation 
volume being probed. While the B/C analysis leaves a 
degeneracy between the diffusion strength and L, nomi- 
nally allowing values as small as L ~ 1 kpc, considerably 
larger values (~4-10 kpc) are preferred when also taking 
into account radioactive isotopes J2f|, gamma rays jioj ]. 
cosmic-ray electrons [!l| or radio data [H, [43| . 

We will therefore mainly use the reference ('KRA') 
model for cosmic ray propagation from a recent compre- 
hensive study 44j, featuring L = 4 kpc. As the most op- 
timistic model in terms of constraining DM annihilation, 
we will also refer to the 'CON' model of the same analysis 
(with L = 10 kpc). We use DarkSUSY to calculate 
the propagation of primary antiprotons and the expected 
flux at the top of the atmosphere (building on the proce- 
dure described in Refs. [3,|4B|). In order to derive limits 
on the annihilation cross-section, we then demand that 
the minimally expected astrophysical background of sec- 
ondary antiprotons (which we take from Ref. j47|) plus 
the DM signal do not overshoot the antiproton measure- 
ments by PAMELA (4^ by more than 3c in any data 
point. For our purpose, the only relevant data points 



turn out to be those between 7 and 26.2 GeV; solar mod- 
ulation and convection affect the antiproton flux only at 
energies smaller than about 1 GeV and the associated 
uncertainties therefore have no impact on our limits. 

For comparison, we show our results in Tab. lIIll for the 
same annihilation channels as for the gamma-ray con- 
straints. In principle, final state radiation of electroweak 
gauge bosons in the case of lepton final states would also 
produce antiprotons [4!| - albeit at rather low rates for 
the small values of m x we are interested in here. We 
thus do not include those final states because the result- 
ing constraints would anyway be much weaker than from 
gamma rays or synchrotron radiation (see below) . 



C. Synchrotron radiation 

Let us now turn to the synchrotron radiation emit- 
ted by electrons from WIMP annihilation at or near the 
galactic center. We use a semi-analytical approach to 
derive constraints from radio observations, in particular 
the 50mJy upper limits on the flux density at 408 MHz 
1501 . and summarize our results in Tabic HTH 

We first note that in a magnetic field of strength B 
the peak contribution to synchrotron emission at a given 
frequency v comes from electrons and positrons with an 
energy roughly twice the critical energy, 



E p (v) ~ 0.29~2£ c 



3 • 0.29 eB 



0.46 



GHz 



B 

mG 



(14) 



GeV. 



As discussed in Ref. [52( , close to the galactic center the 
electrons and positrons are deeply in the diffusive regime 
on time scales over which they loose most of their en- 
ergy due to synchrotron radiation. This becomes clear in 
quantitative terms when comparing the energy loss time 
for synchrotron emission t\ oss and the diffusion coefficient 
D evaluated at the electron energy E p (u), 



D 



3 

4 \7re 

2 / 7T 

3 



1/2 



1 



1 



1/1/2 £3/2 
\ ^1/2 



£3/2 



(15) 
(16) 



Here, as usual, the diffusion coefficient is one third of 
the effective scattering length on the magnetic inhomo- 
geneities which was approximated by the gyration radius, 
corresponding to the Bohm limit. Since D/t\ oss ~ r e v ~ 
4 x 10~ 15 (i//408MHz), where r e = e 2 /m e is the clas- 
sical electron radius, we can use the diffusion approxi- 
mation to estimate the length scale Zdiff over which elec- 
trons propagate during their energy loss time. This gives 
Zdift ~ (D iioss) 1 ^ 2 - To quantify this we have to estimate 
the radial profile of the galactic magnetic field strength. 
Assuming energy equipartition between the plasma and 
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general 
WIMP 


cont. gamma 
(dwarfs) 


cont. gamma 
(GC) 


antiprotons 
('KRA', L = 4kpc) 


antiprotons 
('CON', L = lOkpc) 


synchrotron 
(full cone) 


synchrotron 
(r < 1 pc) 


bb 


7.6 (8.1, 8.6) 


21 (22, 23) 


10.4 (11.6, 11.5) 


4.2 (4.7, 4.7) 


27.5 (29.6, 31.2) 


89 (101, 110) 


T + T - 


16 (18, 20) 


14 (15, 16) 






25.8 (29.6, 32.7) 


369 (441, 500) 




145 (168, 190) 


28 (28, 29) 






18.2 (21.8, 24.7) 


427 (515, 589) 


e + e~ 


89 (104, 118) 


14 (11, 13) 






16.1 (19.4, 22.2) 


419 (506, 579) 


W+W~ 


11 (12, 12) 


24 (24, 26) 


9.3 (9.5, 9.8) 


3.8 (3.9, 4.0) 


29.7 (32.5, 34.7) 


122 (139, 152) 



Table III: This table shows the maximally allowed DM annihilation rate (in units of 10 _26 cm 3 s _1 ) into the stated channels, 
assuming a DM mass consistent with a 77 (yZ, 7/1) interpretation of the line observation. See text for details about the various 
limits and underlying assumptions. 



the magnetic field and constant accretion onto the central 
black hole (BH) mass, one obtains [5l| 

B(r) = 7.2 mG x I (i? acc /r) 2 R acc <r< 100i? acc , 
[l0- 4 r>100i? acc 

(17) 

where i? acc ~ 0.04 pc corresponds to the accretion radius 
of the BH. 

The second radial range of Eq. (fl~7]) . where B oc r~ 2 fol- 
lows from the assumption of magnetic flux conservation, 
turns out to be the most relevant for the synchrotron 
emission considered here. Noting that in the vicinity of 
the galactic center, the dark matter density scale height 
l x = \p x {r) J ' p' x {r)\ equals (r s /2)(r/r s ) 1 ~ a for the Einasto 
profile ([T3)> . the ratio ldis/l x becomes 



^diff 



0.3 I — 



2+a 



(18) 



pairs above the peak energy E p that is produced per an- 
nihilation and is obtained numerically from DarkSUSY 
for the various annihilation channels considered here. 

In Table IIII1 wc report the corresponding bounds on 
av that result from the observed 50mJy upper limit on 
the flux at 408 MHz. We note that the bound obtained 
for the truncated integration region is clearly very con- 
servative since there are also contributions to the total 
synchrotron flux in the line of sight that come from radii 
larger than the cutoff distance 1 pc. These contribu- 
tions, however, require in principle the full solution of 
the diffusion equation and are therefore notably harder 
to compute. Simply integrating Eq. (|19[) over the full 
observation volume, on the other hand, results in some- 
what optimistic constraints as it overestimates the actual 
density of the electrons that produce synchrotron radia- 
tion. Nevertheless, this is the typical approach taken (see 
e.g. Refs. [53|, [HI) and we will adopt it in the following 
in order to allow for a simple comparison. 



For distances r < lpc, electrons and positrons thus ba- 
sically loose their energy in situ which implies that the 
local synchrotron emission power is just proportional to 
the local annihilation rate into electron-positron pairs. In 
contrast, at radii r > lpc diffusion washes out the ra- 
dio emission profile compared to the local pair injection 
distribution. 

Carrying out the same procedure as in [52l [53[ , where 
diffusion is neglected and the monochromatic approxi- 
mation on single-electron synchrotron spectra Eq. (|14l) is 
used, we then obtain the following formula for the total 
synchrotron flux density for a Majorana WIMP (an ad- 
ditional factor of 1/2 is necessary if the WIMP is a Dirac 
fermion) : 



F„ 



1 dW s , 



(av) 



4nR% 



du 



E p p x (r)N e (E p )dV , 
(19) 

where we consider two different integration regions: the 
first is a cone with half-a pert ure of 4", corresponding to 
the 408 MHz observation [50| , while the second is defined 
as the intersection of this cone with a sphere of radius 
1 pc around the GC (which roughly sets the limit where 
Eq. (fT9"|) is strictly valid). N e (E p ) is the number of e^ 1 



IV. IMPLICATIONS FOR MODEL BUILDING 

For an Einasto profile as in Eq. (fl"3"|). the annihilation 
rate to two photons that is required to fit the observed 
line signal is given by 



M 77 = 1.27±°;£ x lO^cmV 1 , 
which translates into (crv)^ z = 3. 14^99 x 10~ 27 c 



(20) 



and (<rv) jh = 3.6Stili x Kr 27 cm 3 s _1 for the case of 
annihilation into 7Z and jH , respectively [ll[ . Assuming 
that the line signal is caused by DM annihilation, we 
can now proceed to combine this information with the 
generic relation between tree- and loop-level annihilation 
rates discussed in Section [TT| and the constraints on the 
tree-level annihilation rate summarized in Table IIIII 

For the case of Majorana DM (see Table HJ, the result- 
ing constraints on how much the imaginary part of the 
amplitude may contribute to the loop process is shown in 
Table ITVl in order to be conservative, we always took the 
lower range of the allowed annihilation rate in Eq. ([20]) 
when deriving these values. As stressed before, these 
constraints assume that the loop signal is dominated by 
only one species of intermediate particles. 
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Majorana 
WIMP 


cont. gamma limit 
(GC) 


antiproton limit 
('KRA', L = 4kpc) 


synchrotron limit 
(full cone) 


expected value 
(effective operator) 


bb 


1.0 x 10" 5 (3.1 x 10" 6 ) 


5.1 x 10~ 6 (1.6 x 10~ 6 ) 


1.4 x 10" 5 (4.2 x 10~ 6 ) 


0.44 (0.16) 


T + T - 


4.9 x 1CT 5 (2.0 x 10~ 7 ) 




9.1 x 10~ 5 (4.0 x 10" 7 ) 


0.33 (0.11) 




8.7 x 10^ 7 (3.3 x 10~ 9 ) 




5.7 x 10~ 7 (2.6 x 10~ 9 ) 


0.15 (0.044) 


e + e~ 


2.9 x 10 -11 (8.7 x 10~ 14 ) 




3.3 x 10" 11 (1.5 x 10~ 13 ) 


0.056 (0.015) 


W+W~ 


0.037 (0.21) 


0.014 (0.083) 


0.046 (0.28) 


0.049 (0.067) 



Table IV: This table shows the maximally allowed relative contribution r of the imaginary part of the amplitude to the total 
annihilation rate of Majorana DM particles into 77 (7Z), as defined in Eq. (|10[1 . if the latter is to provide a viable explanation of 
the observed Fermi gamma-ray line at 130 GeV. The imaginary part is assumed to be dominated by the indicated SM particles 
running in the loop. For comparison, the right column shows the expected value of r when the gray blob in Fig. [T]is replaced by 
an effective operator for a 3-point interaction between the respective pair of SM particles and a pseudoscalar that represents 
the annihilating DM particle pair; note that for the case of yZ final states dominated by fermion loops, this even includes the 
'second cut' in Fig. [I] (which arises because the decay Z —> ff is possible and which is not included in our constraints - see the 
discussion at the end of Section HT]) . 



Obviously, fermionic loops are subject to the tightest 
bounds, especially for small fermion masses m,f. Such a 
behavior should be expected from the helicity-suppresscd 
annihilation of 77 or 7Z into a fermion pair with (L, S) = 
(0,0), which results in cr/^/ oc rnj in Eq. We note 
that this scaling, taking into account also r cx Q 4 , can 
be used to derive very good approximations also to the 
limits on r for other quark pairs qq than shown here; 
this is because the fragmentation functions to photons, 
antiprotons and electrons are very similar to the bb case, 
leading thus to essentially the same constraints on the 
tree-level annihilation rate. Charged gauge bosons do 
not exhibit such a suppression and are therefore not as 
strongly constrained; still, the imaginary part of a W 
loop with 77 (jZ) final states may at most contribute 
~1% (-10%) to the total rate. 

Let us mention in passing that two-gluon final states 
provide very stringent complementary constraints if the 
loop is dominated by quarks [57j • For example, we expect 
(av) gg = 2a^/(9Q 4 aU(av) 11 > 3 x 10- 26 Q- 4 cm 3 s- 1 
in case the observed line corresponds to a 77 final state 
[5^ . Using the fact that both antiproton and continuum 
gamma-ray limits for gg final states are very similar to 
the ones for bb listed in Table HTT1 this implies that down- 
type quarks (Q = — 1/3) in the loop are clearly excluded 
for a 77 signal while up- type quarks (Q = 2/3) may still 
remain marginally consistent with those constraints. If 
the 130 GeV feature is due to a 7Z or jh loop, however, 
the corresponding annihilation rate into gg depends on 
the loop topology and no model-independent conclusions 
are possible. 

For comparison, wc also calculated the full loop ampli- 
tude assuming that the 4-point interaction between DM 
and intermediate particles / in Fig. [T] can be replaced 
by an effective operator for a 3-point interaction between 
I and a pseudoscalar particle that represents the initial 
DM pair (see also Appendix IA 2[) . This corresponds to 
the limit where all new physics states heavier than the 
DM particle are integrated out and therefore cannot con- 
tribute to the loop signal. If, as we always assume, only 



one type of SM particles contributes to the loop, this limit 
does not change the result for the imaginary part of the 
amplitude; it may, however, in principle underestimate 
its real part. Wc report the numerical result of this cal- 
culation (using FeynCalc [55| and LoopTools [5(|) in the 
last column of Table [TVl as can be seen, it it always larger 
than the corresponding upper limits reported in the same 
Table (apart from the case of 7Z final states mediated 
by a W boson loop, which is marginally compatible with 
antiproton limits). 

Even if also heavier particles contribute to the loop, 
which is the more commonly encountered situation in 
typical WIMP scenarios, the loop signal gets strongly 
enhanced if the virtual particles can (almost) be put 
on shell. The large annihilation rates necessary to ex- 
plain the observed signal therefore seem to require rather 
generically that not only the real part, but also the imag- 
inary part of the amplitude provides a significant contri- 
bution - implying that for large annihilation rates, the 
values of r in realistic models should not be too far be- 
low those obtained in the effective operator approach. In 
that sense, our bounds present a serious challenge to any 
model-building attempt trying to attribute the observed 
130 GeV feature to a loop that involves SM particles, with 
the only possible exception of top quarks. 1 

For illustration, let us have a closer look at a particu- 
larly well-motivated example for a Majorana DM candi- 
date, the lightest supersymmetric neutralino [l8j |. which 
is a linear combination of the superpartners of the gauge 
and Higgs fields, 

X = x°i = N n B + N 12 W 3 + N 13 H° + N U H% . (21) 



Top loops obviously do not have an imaginary part because 
mt > m x . However, the contribution from top quarks to the 
loop is often connected to the contributions from other quarks 
in a simple and unique way (e.g. through the Yukawa coupling 
with an s-channcl Higgs); the imaginary part of those loops may 
then be used to constrain even such scenarios. 
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Figure 2: Left: For a large scan over MSSM and cMSSM models with m x ~ 145 GeV, this figure shows the loop-level annihilation 
rate (av)-,z vs. the relative contribution r of the imaginary part of the amplitude to this cross section. Models where the 
neutralino is dominantly a Higgsino (Bino) are indicated by blue (red) symbols; green symbols refer to mixed neutralino DM. 
Filled circles correspond to models where thermal production leads to the correct DM density today, while upper (lower) 
triangles indicate a too large (small) relic density. The shaded area contains only cMSSM models in the co-annihilation region 
(see text for further details). Also shown for comparison are the most relevant limits from Tab. lIVI Right: Same, for neutralino 
annihilation into 77 and m x ~ 130 GeV. 



For this purpose, we plot in Fig. [5] the ratio r as a func- 
tion of the loop annihilation rate av for a large num- 
ber of supersymmetric models that resulted from a scan 
(for details, see Ref. Q) over the parameter space of 
the cMSSM and a phenomenological MSSM-7. Here, 
we keep of course only models with the correct neu- 
tralino masses to account for the observed 130 GeV fea- 
ture. In the figure, we also indicate whether the neu- 
tralino is mostly Bino (Zg = |iVn| > 0.9), Higgsino 
(Zjj = \N 13 \ 2 + \N 14 \ 2 > 0.9) or mixed, and whether 
thermal production leads to the correct relic density. 

As a first remark, one can clearly see that it is es- 
sentially impossible to explain the required large anni- 
hilation rate with thermally produced neutralinos. For 
the case of annihilation into jZ (5t|, furthermore, we es- 
sentially recover our general expectation outlined in the 
preceding paragraphs: in order for the loop-signal to be 
large, there must be a sizable contribution r from the 
imaginary part of the amplitude. This confirms that our 
limits are as stringent as advertised before. The largest 
rates, in particular, are obtained if the neutralino has a 
considerable Higgsino fraction; in this case the annihila- 
tion rate into "/Z is dominated by H^-boson loops and 
we include, for comparison, the corresponding limit from 
Tab. [TV] in Fig. [2] 

For the 77 amplitudes [60], however, there is a large 
number of models with (av) ^ ~ 10 _29 cm 3 s _1 that do 
not follow this general expectation and instead show an 



a) b) 




Figure 3: For small mass differences between sfermions / and 
Bino-like neutralinos \, diagram b) completely dominates the 
process \X ~~ > 77- This diagram, however, only has a real 
part and the largest contribution to the imaginary part of the 
amplitude derives from diagram a) - though W boson loops 
start to contribute with about the same size already for a 
relatively small Higgsino fraction of the neutralino. 

unexpectedly small value of r. While such small anni- 
hilation rates are maybe not too relevant in terms of a 
possible explanation of the Fermi observation, it is nev- 
ertheless quite instructive to discuss the origin of this 
feature in some more detail. The models in question lie 
exclusively in or near the so-called coannihilation region 
of the cMSSM, where the lightest neutralino is an essen- 
tially pure Bino that is almost degenerate in mass with 
light slcptons (in particular the f). For these models, 
slcpton box diagrams (shown in Fig. [3]) dominate the am- 
plitude, with the by far larger contribution coming from 
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diagram b) that receives a significant enhancement for 
small mass differences mj — m x . This diagram, however, 
does not have an imaginary part because the virtual lep- 
tons cannot simultaneously be put on shell. While this 
already explains why r should be considerably smaller 
than expected from our general discussion, there is an- 
other effect which can further decrease its value by up to 
4 orders of magnitude or so: in some cases, the imaginary 
part of the amplitude receives an almost equal contribu- 
tion from slepton box diagrams (i.e. the left diagram in 
Fig. [3]) and W boson box diagrams - albeit with a differ- 
ent sign and thus potentially leading to large accidental 
cancellations. As indicated in the figure, some of the 
models in the co-annihilation region therefore evade our 
bounds on lepton loops even when electrons contribute 
significantly. The largest and thus most relevant annihi- 
lation rates, on the other hand, are again obtained for 
Higgsinos; in this case, the rate is dominated by W bo- 
son loops and thus excluded by our constraints (see also 
Refs. OS El). 

The above example of ncutralino DM nicely illustrates 
the general caveats one should keep in mind for our anal- 
ysis: Much smaller values for r than the limits presented 
in Table IIVI can be realized in concrete models if 

• there are cancellations between several channels 
contributing to the imaginary part of the ampli- 
tude or 

• the imaginary part for the dominant channel is sup- 
pressed with respect to the real part by some sym- 
metry of the initial state that can only be satisfied 
for diagrams without imaginary contributions. In 
the ncutralino case discussed above, e.g., the pro- 
duction of intermediate states in the left diagram 
of Fig. [3] is helicity-suppressed, while this is not the 
case for the right diagram. 

Note that the first point only formally avoids our con- 
straints as there is of course no such cancellation for the 
tree-level final states and thus for the production of sec- 
ondary particles. The second point, on the other hand, 
does present a fundamental limitation to our analysis as 
it corresponds to the situation when the dominant contri- 
bution to the loop signal does not show the same topology 
as displayed in Fig. [TJ 

In this particular case, a sufficiently enhanced line sig- 
nal mediated by SM loop particles could thus in prin- 
ciple explain the observed 130 GeV feature without vio- 
lating constraints from tree-level annihilation products. 
However, it is exactly in this situation that internal 
bremsstrahlung (IB) from charged virtual particles gen- 
erally proceeds at even larger rates. Majorana DM par- 
ticles Xj e.g., annihilate with much larger rates into ff"/ 
final states than into the (helicity-suppressed) // at tree 
level [6l| or the (a om -suppressed) 77/7Z 0], at least for 
small mass differences between \ an d /■ The crucial 
point is that this type of 'virtual' IB [|| gives rise to an 
equally pronounced spectral feature in gamma rays, in- 
distinguishable from a line with the energy resolution of 



Fermi. It is thus not surprising to see that the ratio of 
signal to secondary photons can be much larger for IB 
than for line photons [ll|, which evades the constraints 
discussed in this article and may be argued to strengthen 
the case for a possible IB explanation [H, HH, Gjl, E3 of 
the signal. 



V. CONCLUSIONS 

The recently found evidence for a line-like spectral fea- 
ture at 130 GeV in gamma-ray data towards the galactic 
center 0, Q has triggered an enormous activity both in 
terms of independent data analyses and possible explana- 
tions for the observed excess (see Ref. [Ill] for a review). 
The most exciting and far-reaching explanation would 
be in terms of annihilating DM particles \- I n this ar ~ 
ticlc, we have studied the possibility that the signal is 
produced by monochromatic photons resulting from the 
loop-suppressed process XX ~^ 77 7 1% or 7^j an d consid- 
ered possible consequences for DM model-building. 

By using the optical theorem, in particular, we could 
show that the possibility of standard model particles run- 
ning in those loops is highly constrained by both con- 
tinuum gamma ray, antiproton and synchrotron cosmic 
ray data. In fact, it was noted before that the anni- 
hilation rate required to explain the 130 GeV feature is 
considerably larger than typically expected for WIMPs, 
so that the associated production of standard model par- 
ticles at tree-level would likely be in conflict with both 
continuum gamma-ray and antiproton cosmic-ray obser- 
vations. Here, we have made this argument much more 
rigorous by adopting a systematic treatment in an as 
model-independent way as possible. 

We note that our results provide important constraints 
for any model-building attempt to explain the 130 GeV 
excess. In particular, a direct coupling between DM and 
standard model particles (and new, heavier particles) is 
something very generic in all realistic WIMP frameworks 
that arise in scenarios that simultaneously try to address 
the fine-tuning problems of the standard model. In that 
sense, our findings have far-reaching consequences and 
are by no means restricted to the case of neutralino DM, 
which we have explicitly shown to be essentially impossi- 
ble to reconcile with a monochromatic line interpretation 
of the data. While it is in principle possible that the loop 
signal is dominated exclusively by new particles that are 
heavier than the DM particle, such an option may be 
more difficult to motivate from a global model-building 
point of view. 

DM annihilation via internal bremsstrahlung, on the 
other hand, is a promising alternative explanation: given 
that it may proceed at much larger rates than the tree- 
level annihilation into standard model particles, it is 
presently not constrained by the data. One should also 
note that for neutralino DM, e.g., there is only a factor 
of a few missin g in the annihilation rate to explain the 
observed excess 11| - something which might be possible 
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to achieve by a larger normalization due to a higher local 
DM density than assumed here [(33[ or an increased DM 
density near the galactic center due to the presence of 
baryons, as found in a very recent simulation of a Milky 
Way-like galaxy [(34| • It will thus be exciting to await ex- 
periments with improved energy resolution and/or statis- 
tics^ GAMMA-400 [H or HESS-II [H, that will be 
able to distinguish between the peculiar spectrum of in- 
ternal bremsstrahlung and a monochromatic line (29j . 
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Appendix A: Optical theorem revisited 

In this appendix, we provide for convenience some 
practical tools for the actual calculation of the ratio be- 
tween loop- and tree-level cross section. In particular, we 
describe two rather different strategies which, however, 
of course lead to consistent results. 



J 



1. Standard model rates in Helicity and (L.S) basis 



The first approach is to actually compute the SM cross 
section £>/_>./ as it appears in Eq. ([5]), where / is the par- 
ticle species supposed to dominate the imaginary part of 
the loop and / = 77, jZ, 7/1 is summed over all polariza- 
tion states with the same J p value as the initial state of 
two DM particles (recall that we are only interested in the 
limit of vanishing relative velocity of the DM particles, 
so that J is simply given by their total spin). 

This can be done in different ways, for example by 
expressing (!/_>./ in terms of final-state plane waves as 



§§^7 Yl / d2 P/l- M (/,P/,S / ,m S/ )^(/,L,S)(s)| 2 , 

ms } 

(Al) 

where the quantum numbers (J, M, L, S) of the in- 
termediate state / are fixed. The amplitudes 



M 



(f,p f ,S f , mSf )~>(I,L,S) 



can then be expressed in terms of 



the plane wave amplitudes M( It p t s,m s )^-{f,pf,S f ,m Sf )(s) 
with the aid of Clebsch-Gordan coefficients and spherical 
harmonics, namely 



M 



(f, Pf> S f ,m Sf )^(I,L,S) 



(s) 



= Yl d 2 pY L , mL (p)(J,M\L , llli, S, ms)lM.(f,p f .s f ,ms f )^(I,p,S,ms)( s ) 
m L ,m s 

= Y J d 2 pY L , mL (p)(J 7 M\L,m L7 S 7 ms)i(S,m s \S l7 m Sl ,S 2 ,ms 2 )i 



X - /V/J (/,P/,S , / ,m S/ )^(/,p,Si,rn Sl ,S 2 ,ms 2 )( S ) i 
I 



(A2) 



where we defined the spin and its projection of each in- 
termediate particles as (S^msJ with i = 1,2. 

For J = (realized for both Majorana and Scalar DM), 
on the other hand, it proves useful to consider instead 
the helicity basis {| J, M, Ai, A2)}. In the following, we 
will provide details for this approach. Transformation 
rules between the (L, S) and helicity bases, and a set of 
interesting properties of the latter, can e.g. be found in 
Ref. [68| . Following those conventions, in particular, one 
has 



{J,M,L,S\J,M,X 1 ,X 2 



(A3) 



2L + 1 



J, M = X\L, mj, =0, 5, ms = X) 



2J + 1 

x (S 1 ms = X\Si,ms 1 =X 1 ,S 2l ms 2 = -X2) , 



where A = A x — A 2 and the < ...|... > on the right-hand 
side are the standard Clcbsch-Gordon coefficients. 



a. Determination of final and intermediate states 

J = two-photon states are spanned by just two he- 
licity eigenstates | + 1,+1) and | — 1,— 1), with 

1/5) = -±=| + 1,+1> + -^1-1-1) (A4) 

I/m) = --±=| + 1,+1> + -5=|-1-1) (A5) 

being the CP-even and CP-odd eigenstates, respectively. 
The only allowed two-photon final state |/) for scalar 
(Majorana) DM is thus \fs) (|/m))- The same is true 
for 7Z in the final state, while 7/1 final states are not 
possible for J = (note that the combination L = 1, 
S = 1 cannot be constructed with a massless photon and 
a scalar particle). 

Let us now turn to the possible two-body intermediate 
states /. In case of fcrmions, there are only two indepen- 
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dent hclicity eigcnstates with J = 0, namely | + §,+§) 
and | — 5,— §)) and the CP eigenstates coincide with the 
partial wave basis elements: 



rff\ 



= \L = 0,S=0) 



\L = 1,S=1) 



1,1 1 

V5 l + 2' + 2> 



-I -I) 
2' 2 ; ' 

(A6) 

I -I) 
2' 2 ; ' 

(A7) 



In the case of two W bosons, on the other hand, there are 
three independent states with J = 0; one CP-odd state 
(given by \L = 1, S = 1)) and two CP-even states which 
can e.g. be chosen as {\L = 0, S = 0), \L = 2, S = 2)} or, 
cquivalently, transverse and longitudinal states: 



t wwm 
L S / 

T ww.u 



I T WW\ 
\ 1 M I 



|0,0), 

\L = 1,S=1) 



-l.+l) 



(A8) 
(A9) 

1,-1). 

(A10) 



Which linear combination of 



WW,t\ 



and 



I T WW,l\ ■ 

\I S ' ) is re- 



alized as the dominant contribution to the loop process 
(in the case of scalar DM) can thus not be determined 
in terms of the symmetries but depends on the SM and 
DM interaction. 



b. Computation of helicity amplitudes 

With the above definitions, it is now straight-forward 
to express any amplitude (I\T\f) with J = as a sum 
of helicity amplitudes M\ . Xi = (Aj, Ai|T|A 7 , A 7 ) = 



-Xj- 



Those M\ 



can also be written in terms 



of plane waves with spins pointing in the direction of 
motion, i.e. the familiar 4-momcntum two particle states 
with definite helicities. For the states, we have [63 



\J,M, Ai,A 2 ) 



(All) 



2J+ 1 

47T 



d 2 kd J M ,{ 



where again A = Ai — A2, and (8', if') are the angles defin- 
ing fe with respect to some arbitrary cartesian system of 
coordinates. By noting that = 1, we may thus write 



d 2 fc 



(J = 0,M = 0,A/, A/|T|fe; A 7 ,A 7 

4vr(J = 0,1/ = 0, A/,A 7 |T|£; A 7 , A 7 ) 
d 2 p(p; A / ,A / |T|£;A 7 ,A 7 ) 



= 2tt 



Mx.. Xi {0) sm6 AO , 



(A12) 



where the second line follows from the spherical sym- 
metry of the state | J = 0, M = 0, A/, A/) and rotational 
invariance of the interaction, and in the last step we used 
the fact that the integrand does not depend on the az- 
imuthal angle ip (z is an arbitrary direction, taken to be 
aligned with the photon momentum, and 9 is the interac- 
tion angle with respect to it, i.e. the center-of-mass angle 
between the different momenta/hclicitics). 

■M-\ 1 -\ I {9), as introduced above, is the amplitude for 
two photons (or one photon and one Z boson) with equal 
helicities and opposite momenta that annihilate into a 
pair of charged particles with opposite spins parallel to 
their momenta. These amplitudes can directly be cal- 
culated by standard Feynman rules (see below for our 
conventions and more details) and we find 



KA n^w + w- im _ 2e 2 (2A 2 y + 2AtyA 7 /3/ - 7f 2 ) 



= jZowl-ftcaee ty [ (1 " im Z ,2m x f^e )Pl + 2/3/ A 1/2 A 7 ] + 

+ A 7 C^(A + 2A 1/2 A 7 )/3/cos0) , (A15) 

o Q ^/(A 7 + X w Pi) 2 cos 2 9 + (A 7 /3/ + X w Pi) 2 sin 2 9 
Ml M0 (9)=2egcos9 w M l - ft cos 2 9) ' (M6) 

U 7Z^W+w-,^ 0t , nrna a PA 1 - m z/ 2m w) cos2 g + (1 ~ ml/Ami, - m|ffi/8m^) sin 2 

^A 7! A W =0 (e)—2egco*e w , (A17) 



r 



where the reduced velocity /3j of the intermediate state particles I is defined in Eq. ((5]), with a corresponding 
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definition holding for the reduced velocity /?/ of the fi- 

1/2 

nal state particles, and ji = (l — 0f) = m x/ m f(w)- 
Furthermore, we introduced = T 3 - 2Q sin6*n/ and 
= T3, where Q and T3 are the fermion charge and 
weak isospin projection (for left-handed fermions), re- 
spectively. From these expressions, the helicity ampli- 
tudes At£ are obtained by integration over 6. As 
explained above, the amplitudes (I\T\f) that appear in 
Eq. can then easily be computed as a sum over all con- 
tributing M\ . Xi , leading to the results shown in Tab.|U 



c. Conventions for polarization vectors and spinors 



In arriving at the amplitudes (|A13IIA17|) . we adopted 
the following explicit representations of the appearing 4- 
momcnta (with p/E p = /?/): 

fc 7 = (k,0,0,k) , (A18) 

k %z) = {E k ,0,0,-k) , (A19) 

Vi = {E p ,psm6,O,pcos0,) , (A20) 

pj = (E p ,-psm9,0,-pcos9) . (A21) 

For fermion states J, we used the Weyl representation for 
gamma matrices and the following set of spinors: 

/ ^/E p - 2Ai /2 pcosf 4 ' ) 

^(2B + (l-2\ 1/2 )-n s 



u(Ai/ 2 



yJE p - 2A 1/2 psin(- 
y/Ep + 2A 1/2 pcos( 4 j 

TTv • / 20+(l-2A 1/2 W\ 

y^/E p + 2\ 1/2 psml K 4 J J 



(A22) 



«(-Ai/ 2 ) 



, 2e-(l-2Ai/ 2 )7r\ \ 

v / J B p -2A 1/2 psm( 4 u '' ' \ 



' 29-(1-2A 1/2 )tt s 
> 26>+(3+2A 1/2 )ir < 



^/E p - 2Ai/ 2 pC0S 

^/Sp + 2Ai/ 2 psin 

/-S, — rrr /20+(3+2A 1/2 )ir" 

^ V% + 2A 1/2 pcos( 



.(A23) 



/ 



Last but not least, the polarization vectors for the vector 
bosons in our conventions read 



e 7 (A 7 ) = rO,-A 7 /v / 2,-i/v / 2,o) 
^(Z)(\) = (0,A 7 /V2,-i/V2,0^ 



\Y ) 



(1 - A^) 7 //?i 



JL + (1-X 2 w hi sini 

A W ^2 



72 



(A24) 
(A25) 

(A26) 



V ^f^ + (i-A 2 M/ ) 7 /cos0 y 



-(1 - A^) 7jj 8j 



+ (l-A 2 y ) 7/ sin6' 



A]£ COS £ 

V2 

\ 2 i 
A ^V2 

_A 2 ^ if + (1 _ A 2 /)7/COS0 y 
2. Effective coupling approach 



V 



(A27) 



The second option (used e.g. in Ref. [28|]) is to intro- 
duce a fictitious particle 0, with = 2m x and the 



same parity and spin as the initial state. One then needs 
to specify an effective 3-point coupling between <f> and 
the two intermediate states /, which corresponds to con- 
tracting the 'blob' in Fig. [T] to a point and considering 
the decay of <p instead of the annihilation of two particles 
X- While this approach does not explicitly make refer- 
ence to the general formulae presented in Section HH it is 
sometimes more convenient to follow in practice. Note, 
in particular, that our general discussion in Section |TT] 
implies that the ratio r^r between loop- and tree-level 
cross section, or decay rate, does not depend on the form 
of the effective operator that is chosen (as long as it cor- 
rectly connects particles of the required spin and parity). 

The reason why this approach can be computationally 
simpler is that one can directly compute the processes 
i — > I and i — > / by standard means, i.e. by using the 
full polarization and/or spin sums and without having to 
worry about explicit projections on specific (L, S) states. 
Even the loop integral for i — > / contains only three 
legs and is thus relatively straight-forward to compute 
by standard techniques. In fact, one may even directly 
compute Disc M.i^j = 2iQvWi_>./ by using the cutting 
rules proven by Cutkosky Jg^: simply take the full ex- 
pression for M.i^f and replace 



1 



-2irid(q 2 - m 2 j) 



(A28) 



for the two propagators of the (necessarily on-shell) inter- 
mediate particles. This 'trick' reduces the loop integral to 
two angular integrals (one of which is trivial) and makes 
the integrand much simpler because only one remaining 
propagator is involved. Note that the effective operator 
approach also allows to rather easily calculate the con- 
tribution to the imaginary part of the loop integral that 
derives from cutting along the dotted line of Fig. [T] which 
appears only for 7Z final states dominated by fermion 
loops and which is neglected in Eq. ^ and thus in the 
approach described in the previous subsection. 

For completeness, let us finally list all effective opera- 
tors that are needed to compute the relevant rates for an- 
nihilating Majorana DM. In this case, </> is a pseudoscalar 
that interacts with fermions via 



47/ = 94>J^h 5 f 
and with gauge bosons via 



/»int 

'-4>ww 



= g+,w4>(d li w+){d p w- 



(A29) 



(A31) 



where F£ v denotes the SU{2) field strength. As stressed 
above, the coupling constants g^,... that appear in these 
expressions drop out when calculating ratios like r; alter- 
natively, they may be fixed by normalizing to either the 
tree or loop level rate. 
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